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ABSTRACT 

The stellar content of the open cluster NGC 2437 (Messier 46) is investigated 
using moderately deep u*,g', and r' MegaCam images. When compared with 
solar metallicity isochrones, the (g', v! — g') and (r', g' — r') CMDs are consistent 
with an age log(t yr ) = 8.35 ± 0.15, a distance modulus /i = 11.05 ± 0.05, and 
a color excess E(B — V) = 0.115 ± 0.035. The r' luminosity function (LF) of 
main sequence stars in the magnitude range r' < 17 (i.e. masses > 0.8 M Q ) 
has a shape that follows solar neighborhood star counts. However, at fainter 
magnitudes the cluster LF is flat, in contrast with what would be expected from 
solar neighborhood counts. The clustering properties of stars in NGC 2437 are 
investigated by examining the two-point angular correlation functions of main 
sequence stars in different brightness ranges. Main sequence stars fainter than 
r = 17 are less centrally concentrated than brighter stars and are found over a 
larger area of the sky, suggesting that there is a corona of faint main sequence 
stars around NGC 2437. Based on the flat LF and extended spatial distribution 
of faint stars, it is concluded that NGC 2437 is actively shedding stars with 
masses < 0.8 M . 

Subject headings: Star Clusters and Associations 



1 Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Sciences de l'Univers of the Centre National de 
la Recherche Scientifique (CNRS) of France, and the University of Hawaii. 
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1. INTRODUCTION 

Star clusters are basic targets for observational studies of stellar and galactic evolution. 
However, the majority of clusters are short-lived, being disrupted after only a few Myr (e.g. 
Lada & Lada 2003) probably due to large-scale mass loss driven by supernovae and stellar 
winds. The clusters that survive their infancy are likely those with initial conditions that 
(1) allow them to develop regions early-on where stars are the dominant mass component, 
thereby reducing the disruptive influence of gas loss, and/or (2) are in a contracting phase 
of their dynamical evolution at the time of peak gas expulsion (e.g. Smith et al. 2011). 
Clusters that are seen at the present day may then be biased examples of the majority of 
clusters that initially formed. 

The spatial distribution of stars within clusters provide insights into their past history. 
Clusters that are younger than a few tens of Myr can be surrounded by diffuse halos that 
contain some of the most massive cluster members (e.g. Davidge 2012). The spatially 
dispersed nature of these massive stars is likely the consequence of the large-scale loss of gas 
early in the evolution of the cluster (Bastian & Goodwin 2006). The central bodies of older 
clusters are surrounded by coronae populated by low mass stars (e.g. Kholopov 1969), which 
are expected to move to the cluster periphery due to dynamical relaxation (e.g. Prisinzano 
et al. 2003; Sharma et al. 2006). Stars in the outermost cluster regions will be the most 
susceptible to stripping, and as clusters in the disk disperse they may leave debris trails that 
extend over large swaths of the sky (e.g. Chumak & Rastorguev 2006; Chumak et al. 2009). 

In the present study, u*,g', and r' images obtained with the Canada- France-Hawaii 
Telescope (CFHT) MegaCam are used to examine the stellar content of the open cluster 
NGC 2437. NGC 2437 is an interesting target for investigation as it is nearby and has a 
mass and age that make it pertinent for studies of cluster disruption. Basic properties of 
NGC 2437, taken from the WEDBA database (Mermilliod 1995), are listed in Table 1. 
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It is somewhat surprising that - despite being a relatively nearby and well-populated 
intermediate age cluster - NGC 2437 has not yet been extensively studied with electronic 
detectors. Sharma et al. (2006) use CCD images to investigate stars in NGC 2437 to as 
faint as V ~ 18 — 19, and find a distance of 1.51 kpc (/x = 10.9) and an age \og(t yr ) = 8.4. 
Their data trace the radial distribution of stars out to 20 arcmin from the cluster center. 
Sharma et al. (2006) find that NGC 2437 is more spatially extended in the near-infrared than 
at visible wavelengths, hinting that the redder (lower mass) cluster members may be more 
widely distributed than the bluer (higher mass) members. Still, despite this hint that mass 
segregation may be underway, Sharma et al. (2008) find that stars as faint as V ~ 18 — 19 in 
NGC 2437 have a mass function exponent that is consistent with the Salpeter (1955) value. 

The specific goals of the present study are (1) to estimate the age, reddening, and 
distance of NGC 2437, and (2) to search for signatures of cluster dissruption by examining 
the cluster luminosity function (LF) and the spatial distribution of cluster stars spanning 
a range of masses. A discussion of the observations and their reduction is provided in §2. 
Particulars of the photometric measurements are discussed in §3, while the cluster CMDs and 
LFs are examined in §§4 and 5. The spatial distribution of stars in the cluster is investigated 
in §6. A summary and discussion of the results follows in §7. The use of different criteria 
for identifying objects with suspect photometry and the possible relationship between NGC 
2437 and the planetary nebula (PN) NGC 2438 are investigated in the Appendix. 

2. OBSERVATIONS 

The data were recorded with the CFHT MegaCam (Boulade et al. 2003) as part of 
program 2010BD89. The detector in MegaCam is a mosaic of thirty-six 2048 x 4612 E2V 
CCDs arrayed in a 4 x 9 format with 0.185 arcsec pixel -1 sampling. There are 80 arcsec 
gaps between the three CCD banks, and 13 arsec gaps between individual CCDs. A single 
exposure images a ~ 1 degree 2 field. 
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Details of the observations are summarized in Table 2. Images with short and long 
integration times were recorded to extend the magnitude range of detected stars. The min- 
imum exposure time in u* was longer than in the other filters to avoid potential charge 
transfer problems that could arise from the low sky levels recorded through this filter during 
dark time. 

Initial processing was done with the CFHT ELIXIR pipeline, and this included bias 
subtraction and flat-fielding. The ELIXIR-processed images were corrected for positional 
offsets between images that are recorded in different filters. These offsets vary in size with 
location in the science field, and are typically a few pixels (i.e. < 1 arcsec) near the edge of 
the MegaCam field when u* and g' images are aligned near the center of the detector mosaic. 
Offsets between the g' and r' images are much smaller. These distortions were corrected in 
a differential manner by using the IRAF GEOMAP and GEOTRANS routines to map the 
u* and r' images of each cluster into the reference frame defined by the g' image. 

3. PHOTOMETRIC MEASUREMENTS 
3.1. Stellar Brightnesses 

Target lists and point spread functions (PSFs) were obtained by running routines in the 
DAOPHOT (Stetson 1987) package. There are many unsaturated stars with high S/N ratios 
in the images, and PSFs were constructed by combining 50 - 100 of the most isolated of 
these in each filter. Faint companions were removed iteratively by subtracting progressively 
improved versions of the PSF. 

Stellar brightnesses were measured with the PSF - fitting routine ALLSTAR (Stetson & 
Harris 1988). Sources that are not well-matched by the PSF, and/or that hinder convergence 
to a stable solution are rejected by ALLSTAR. Still, the photometry catalogue that is output 
by ALLSTAR contains sources with poor photometry, and some of these can be identified 
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using information that is provided by ALLSTAR. Following Davidge (2010), such objects 
in the CFHT images were identified using e, the error in magnitude that is computed by 
ALLSTAR. 

Stars that are neither saturated nor blended define a monotonic relation between e 
and magnitude. Objects that depart from this relation tend to be background galaxies 
(e.g. Figure 2 of Davidge 2010), unresolved blends, saturated stars, or cosmetic defects, and 
these were removed from the photometric catalogues. Objects with e > 0.3 magnitudes, 
the majority of which are near the faint limit of the photometry, were also deleted. The 
use of other ALLSTAR-computed indices to identify objects with suspect photometry is 
investigated in the Appendix. 

3.2. Calibration 

Instrumental magnitudes were transformed into the SDSS system by applying linear 
transformation relations with zeropoints obtained from observations of photometric stan- 
dards that are recorded during each MegaCam observing period. The u* — > u' transfor- 
mation warrants some discussion, as the effective wavelengths of the u* and v! filters differ 
by ~ 200 A, in the sense that u* is redder. This difference in effective wavelength is par- 
ticularly acute as the u' filter samples the Balmer discontinuity. The difference between 
the filter passband then affects the ability to obtain reliable v! magnitudes for intermediate 
spectral-type high surface gravity stars from u* magnitudes (Clem et al. 2008). While pre- 
vious experience demonstrates that a reliable u* — > u' transformation can be acheived using 
a linear transformation equation, that result is restricted to intrinsically luminous stars with 
weak Balmer breaks (e.g. Davidge & Puzia 2011; Davidge et al. 2012). 

Bright main sequence stars in NGC 2437 have an intermediate spectral- type, and so 
have strong Balmer discontinuities. Because cluster members have a (more-or-less) common, 
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measureable reddening, their u* magnitudes could be transformed into v! magnitudes using 
- say - a non-linear relation like that found by Clem et al. (2008), which would result in 
corrections of up to ~ 0.4 magnitudes being applied to the u* magnitudes of late B and early 
A-type main sequence stars. However, the application of a non-linear relation that involves 
large corrections to field stars, which have a range of reddenings and evolutionary states, 
is much less secure, as there is the potential to produce final magnitudes that are greatly 
in error. For example, the large transformation corrections from Clem et al. (2008) that 
are appropriate for A and B main sequence stars would be applied in error to stars of later 
spectral-type that have had their reddening underestimated. To prevent introducing such 
large errors into the transformed magnitudes, a psuedo v! magnitude, u'^ was computed for 
each source using a linear transformation relation. In §4 it is demonstrated that the SEDs 
of stars with u'^ magnitudes match those computed using u' over most of the color range 
sampled in NGC 2437. 



3.3. Artificial Star Experiments 

Sample completeness and uncertainties in the photometry - both random and system- 
atic - were assessed by running artificial star experiments. The artificial stars were assigned 
colors and brightnesses that are representative of objects in NGC 2437. As with the actual 
observations, an artificial star was considered to be detected only if it was recovered in at 
least two filters (either u'+g' or g'+r'). The photometric catalogues produced by ALLSTAR 
as part of these experiments were filtered with the e-based criteria that were applied to the 
science data. 

Errors in the photometry were estimated by comparing the input and measured bright- 
nesses of artificial stars, and then computing the mean difference and the dispersion about 
the mean at various magnitudes. Completeness and the uncertainties in the photometry are 
coupled, in the sense that random and systematic uncertainties in the photometry climb 



-7- 



rapidly towards fainter magnitudes when the completeness is below 50%. The magnitude at 
which 50% completeness occurs, which is near v! — g' — r' — 23 in the deep exposures of 
NGC 2437, is thus one measure of the photometric faint limit. 

4. RESULTS: CMDs 

The ((/, u'* — g') and (r', g' — r') CMDs of sources in and around NGC 2437 are shown 
in Figures 1 and 2. Measurements with g' < 14.5 and r' < 14.5 were obtained from the 
short exposure MegaCam images. The open squares are stars with V < 10.5 from Lynga 
(1960), and the published magnitudes were transformed into the SDSS system using the 
relations from Smith et al. (2002) and Jordi et al. (2006). The envelope of ±2cr random 
errors computed from the artificial star experiments is shown in the lower left hand corner 
of each CMD. The uncertainties in the photometry become substantial for objects fainter 
than g' = r' = 20. 

The main sequence of NGC 2437 cuts diagonally across both CMDs, although the 
ability to trace it is compromised at some points by contamination from non-cluster objects. 
These contaminants are a mix of field stars and galaxies, although the contribution from 
galaxies should be modest given the low Galactic latitude of NGC 2437. In addition to 
producing the diffuse spray of points that is scattered over both CMDs, field stars populate 
two distinct features, the most prominent of which is a near- vertical blue plume. This feature 
intersects the main sequence of NGC 2437 near g' m r' ~ 15, complicating efforts to trace 
the cluster sequence in this part of the CMD. Field main sequence stars also populate the 
vertical plume near the faint end of the (r',g' — r') CMD. This feature intersects the faint 
end of the cluster main sequence. 

The (it 1 * — g',g' — r') two-color diagram (TCD) is shown in Figure 3. The cyan line 
is the locus of main sequence stars traced from the — g') and (g',g' — r') CMDs of 
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NGC 2437. The red line is the locus of Galactic main sequence stars from Johnson (1966), 
transformed into the SDSS system using relations from Smith et al. (2002) and Jordi et al. 
(2006), and then reddened to match the E(B — V) of NGC 2437 using the relations in Table 
6 of Schlegel et al. (1998). The agreement between these two sequences suggests that the 
u'^ magnitudes are reasonable proxies for SDSS u' magnitudes over much of the color range 
sampled by these data. 

The spectral-energy distributions (SEDs) of many non-cluster objects will differ from 
those of cluster sources, raising the prospect that at least some contaminating objects can be 
identified based on their location on the TCD. Consider the ridgeline of the field star plume, 
which is shown as the cyan line in Figure 3. This is an almost vertical feature on the CMDs, 
and so occupies a comparatively modest range of colors in the TCD when compared with 
solar neighborhood and NGC 2437 main sequence stars. Still, field stars on this sequence 
are offset ~ 0.1 magnitude along the (g' — r') axis from the cluster sequence on the TCD - 
at least some of the sources on this very prominent feature can be identified based on their 
locations on the TCD. 

For the current study, two samples of objects are defined based on their location on 
the TCD. 'Cluster' stars are those that lie within ±0.275 magnitude in u[ — g' of the NGC 
2437 main sequence on the TCD, while 'non-cluster' objects are those that occupy the rest 
of the TCD. Experiments indicated that this color envelope is the narrowest that preserves 
cluster main sequence stars. The application of a tighter color envelope about the NGC 
2437 sequence on the TCD will reject more non-cluster objects, but at the expense of also 
rejecting cluster main sequence stars. 

The CMDs of objects in the cluster and non-cluster samples are compared in Figures 
4 and 5, where the data have been sorted into three radial intervals that sample equal areas 
on the sky to facilitate an assessment of residual contamination from non-cluster objects. 
While there is residual non-cluster contamination in the cluster sample, the contamination 
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in the upper rows of Figures 4 and 5 is substantially smaller than in Figures 1 and 2. There 
is no obvious binary sequence above the cluster main sequence, suggesting that the fraction 
of binaries with mass ratio near unity is modest in NGC 2437. 

It is evident from Figures 4 and 5 that all but the bluest field plume stars have been 
identified as non-cluster objects. The majority of stars that occupy the plume with r' > 18 
and g' — r' ~ 1.5 in the (r', g' — r') CMD are also identified as being non-cluster in origin. The 
number of non-cluster objects in each radial interval is more-or-less constant, as expected if 
these objects are uniformly distributed on the sky. 

Isochrones with Z = 0.019 from Marigo et al. (2008) are compared with the cleaned 
inner annulus CMDs in Figure 6. The distance modulus and reddening listed in the WEDBA 
database as of mid-2012 have been assumed for this initial comparison. As discussed in §3, the 
MegaCam data are not in the SDSS v! system. Therefore, the isochrones were transformed 
into the v!^ system by finding the offset in u' — g' between the observed and model main 
sequences at a given g' — r' color on the (V — g',g' — r') TCD. This offset was then applied 
to the model v! magnitudes. 

The isochrones in Figure 6 fall above the cluster main sequence, indicating that the 
adopted distance modulus is too small. The distance modulus and reddening were adjusted to 
obtain better agreement with the observations, and the results for fi = 11.05 and E(B— V) = 
0.115, which are adopted for the remainder of the paper, are shown in Figure 7. Comparisons 
with a grid of fi and E(B— V) values indicated that variations in these quantities that exceed 
±0.05 magnitude (fio) and ±0.035 magnitude {E(B — V)) could be identified, and so these 
are adopted as the estimated uncertainties. The reader should keep in mind that these 
uncertainties assume that the stars in NGC 2437 have a known metallicity (Z = 0.019). 

The models shown in Figure 7 suggest an age \og{t yr ) = 8.5 based on the (g f , u'^ — g') 
CMD, and log(V) = 8.2 based on the (r',g' - r') CMD. Therefore, log(V) = 8.35 ± 0.15 
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is adopted for NGC 2437, and additional support for this comes from the bright end of the 
cluster LF (§5). Isochrones with this age do not pass through the two bright blue objects 
observed by Lynga (1960). To the extent that the models are correct, then this suggests 
that these stars (1) are not cluster members, (2) are variable stars, (3) are binaries, and/or 
(4) are not experiencing 'normal' evolution of the type assumed by the models. 

5. RESULTS: LFs 

The r' LF of objects that (1) have g' — r' within ±0.15 magnitudes of the NGC 2437 
main sequence on the (r',g' — r') CMD, and (2) are within 19.3 arcmin of the center of 
NGC 2437 - where the density of cluster stars is highest (Figures 4 and 5) - is shown in 
Figure 8. The faint limit of the LF is defined by the g' measurements, and the artificial 
star experiments indicate that the data are 100% complete for r' < 20, which corresponds 
to > 0.5 M . The LF in Figure 8 has been corrected for non-cluster objects using source 
counts in the outer annulus. 

The r' LF can be divided into three magnitude regimes. At the bright end (r' < 13), 
the LF climbs steadily towards fainter magnitudes. There is then a break in the LF at 
r' ~ 13, which is due to an inflexion in the mass-magnitude relation near ~ 2 M . Finally, 
the LF is more-or-less flat when r' > 16 (i.e. for masses < 1 M Q ). It is demonstrated in 
the Appendix that the flat nature of the LF at these magnitudes is a robust result that is 
reproduced using different criteria to cull sources with poor photometric measurements. 

Model LFs, generated from the Marigo et al. (2008) isochrones and downloaded from 
the Padova database of evolutionary tracks and isochrones web site, |^| are also shown in 
Figure 8. The models assume Z = 0.019, and have ages t = 200 Myr (dashed line) and 
300 Myr (dotted line). A Kroupa (1998) Initial Mass Function (IMF), which is based on 



2 http://stev.oapd.inaf.it/cgi-bin/ cmd 
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solar neighborhood star counts, was adopted. The models have been scaled to match the LF 
between r' = 13.5 and 15.5. 

The change in the character of the LF near r' ~ 13 provides a loose check on the 
distance modulus, and the agreement with the models at this magnitude suggests that the 
adopted distance modulus is correct to within ±0.3 magnitude. While the 300 Myr model 
provides a better match to the observations at the bright end than the 200 Myr model, 
there is a tendency for both models to fall above the NGC 2437 LF at the bright end. Near 
r' = 12.5 the oberved number counts differ from the 300 Myr models at more than the 2a 
level. Better agreement with the model LFs would result with an apparent distance modulus 
that is ~ 0.1 dex greater than deduced from the CMDs, although this would come at the 
expense of degrading the agreement between the isochrones and the cluster main sequence 
on the CMDs at intermediate brightnesses. 

The models depart significantly from the observed LF when r' > 17, with the size of 
the disagreement increasing steadily towards fainter r' . At r' = 20 (main sequence mass 
~ 0.5 M Q ), the number of stars predicted by the model is roughly an order of magnitude 
greater than observed. The LF of the inner regions of NGC 2437 thus differs significantly 
from that of the solar neighborhood when r' > 17 (i.e. < 0.8 M Q ). Possible causes of this 
difference are examined in §7. 

6. CLUSTER STRUCTURE 

The spatial distribution of stars in a cluster contains clues about its past history. The 
spatial distribution and clustering properties of stars in and around NGC 2437 are examined 
in this section using two approaches. The two-point angular correlation function (TPCF) of 
cluster stars is examined in §6.1, while the physical distribution of cluster stars on the sky 
is considered in §6.2. 
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As with the analysis of the LF in §5, objects within ±0.15 magnitude in g' — r' about 
the main sequence ridgeline are adopted as cluster members, although the result of applying 
a narrower extraction window is also investigated. Cluster stars are divided into various 
brightness regimes using the comparison between the observed LF and models in Figure 
8 as a guide. The BMS (bright main sequence) sample includes the most massive main 
sequence stars in the cluster, which have r' < 13. The IMS (intermediate main sequence) 
sample includes objects from r' = 13 to r' = 17. The LF increases slowly towards fainter 
magnitudes in this magnitude interval, at a rate that is more-or-less consistent with that 
expected from the models shown in Figure 8. The LMS (lower main sequence) sample runs 
from r' — 17 to r' — 20. This is the magnitude interval where significant departures from 
the model LFs occur. Finally, a field star (FS) sample is defined that consists of all objects 
that are outside of the ±0.15 magnitude main sequence extraction region. The objects in 
the FS sample are uniformly distributed across the field, and are used to remove signatures 
in the TPCF that arise from detector geometry - such as gaps between detectors and the 
finite size of the MegaCam science field. 



6.1. The TPCF 

The TPCF multiplexes information by considering all possible object-object pairings in 
a sample, and so makes efficient use of the available information. However, as the correlation 
properties of all objects are considered in concert then there is not a direct relation between 
features in the TPCF and location within the area examined. Still, some generalizations can 
be drawn about the locations of groups that contribute signal to the TPCF. Basic geometric 
arguments imply that signal at the largest possible separations must originate near the edges 
of the sampled area. The potential for ambiguity grows as progressively smaller separations 
are considered, as these separations can be measured over a larger fraction of the field. 
Another mitigating factor is the density distribution of stars in a cluster. 
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The distributions of separations obtained for all possible object pairings (hereafter 
'separation functions' - SFs) were calculated for sources in each of the samples. TPCFs 
were then constructed by dividing the SFs of the BMS, IMS, and FMS samples by the SF 
of the FS sample, and scaling the results according to the number of pairings. The resulting 
TPCFs are compared in the top panel of Figure 9. 

The objects in the three samples have different clustering properties. The TPCF of 
the BMS sample has by far the largest amplitude, with power exceeding zero detected out 
to separations of ~ 60 arcmin. An inflexion point in the TPCF of the BMS sample occurs 
near separations of 20 - 25 arcmin (~ 9 — 12 parcsecs), hinting at a possible change in the 
structural characteristics of the cluster 10 - 12.5 arcmin from its center. A discontinuity is 
also seen in the TPCFs of IMS and FMS stars at this separation. 

A complication in interpreting the TPCFs is that field stars are located over a wide 
range of distances along the line of sight, and so binaries with a range of different physical 
separations are included in each angular interval. Consider very wide stellar binaries, which 
have spatial separations ~ 0.04 parcsecs (e.g. Larson 1995; Simon 1997). Binaries in the 
field with separations of 0.04 parsecs would be separated by a few arcmin if viewed at a 
distance of 100 - 200 parcsecs, but similar binaries at the distance of NGC 2437 will have 
separations < 5 arcsec. It is thus of interest to consider the ratios of SFs of objects that are 
in the cluster. 

The ratio of the IMS and FMS SFs is shown in the lower panel of Figure 9. The IMS 
and FMS samples have similar clustering properties at separations < 25 arcmin. However, 
at larger separations IMS stars are less clustered. This is consistent with objects in the FMS 
sample contributing a larger fraction of the objects in the outermost regions of NGC 2437 
than stars in the IMS sample. 

The TPCFs of the IMS and FMS samples have substantially smaller amplitudes than 
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the BMS TPCF. This reduction in amplitude is due - at least in part - to differences in 
contamination from field stars, which are uniformly distributed across the field and cause 
the TPCFs of field star-contaminated samples to be flatter than those of samples that are 
free of such contamination. Star counts in the inner and outer annuli (§5) suggest that as 
many as < 60% of the stars in the IMS and FMS samples may be field objects, while only 
< 10% of the BMS sample may be field stars. 

The extent to which non-cluster objects affect the TPCFs can be assessed by comparing 
the TPCFs of samples constructed using different extraction widths along the color axis on 
the CMD. The fractional contamination from field stars is reduced as progressively narrower 
extraction windows are applied, although this reduction is obtained at the expense of also 
culling cluster stars that scatter about the main sequence on the CMD. To investigate the 
impact of non-cluster contamination on these data, a second set of TPCFs were computed 
for the IMS and FMS samples using a ±0.05 magnitude wide band in g' — r' . While not 
shown in Figure 9, the TPCFs of the IMS and FMS samples constructed with this narrower 
window have amplitudes that are ~ 2x larger than those obtained with a ±0.15 magnitude 
cut. This is consistent with reduced levels of non-cluster contamination. 

The ratio of the SFs constructed with the ±0.05 mag extraction window is shown as 
a dashed line in the lower panel of Figure 9. The ratio of the IMS and FMS SFs obtained 
with the ±0.05 mag window differs from that obtained with the ±0.15 mag window, in the 
sense that the IMS sample is more centrally concentrated than the FMS sample when the 
narrower extraction window is used. This comparison reinforces the result discussed above 
that the cluster stars in the FMS sample are more diffusely distributed than those in the 
IMS sample. 
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6.2. The Projected Distribution of Stars 

Additional insights into the structure of NGC 2437 can be obtained by visually ex- 
amining the spatial distribution of objects. The distributions of sources in the BMS, IMS, 
and FMS samples are shown in Figure 10. The top row shows the source distributions of 
objects that fall within a ±0.15 mag extraction window about the main sequence locus on 
the (r', g' — r') CMD, while the bottom rows shows the distribution of sources within a ±0.05 
mag extraction window. The pixel intensities reflect the number of stars within 1.5 x 1.5 
arcmin gathers. This coarse bin size was adopted to suppress stochastic fluctuations in the 
star counts, still, the BMS distributions are dominated by single objects outside the main 
body of the cluster. 

Features that are present in the top row of Figure 10 also tend to be present in the 
bottom row. That the spatial distributions of each sample type are not sensitive to the 
width of the extraction window on the CMD suggests that prominent structures in Figure 
10 are dominated by cluster stars, and not field stars. The objects in the various samples 
have different spatial distributions, with the BMS and IMS samples being more compactly 
distributed than the FMS sample. This is consistent with the comparisons between the 
TPCFs of the samples in the previous section. 

The IMS and FMS distributions have irregular morphologies at large radii, and the FMS 
distributions show a concentration of objects to the south and west of the cluster center. 
This feature is likely not due to field stars, as it is equally strong in the top and bottom 
rows of Figure 10, indicating that it is populated by objects with photometric properties 
that tightly follow the cluster main sequence. The nature of this feature is discussed further 
in the next Section. 
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7. SUMMARY AND DISCUSSION 

The stellar content of the open cluster NGC 2437 (Messier 46) has been investigated. 
Despite being a moderately massive nearby intermediate age cluster, NGC 2437 has not yet 
been extensively surveyed with modern detectors; the one exception is the study by Sharma 
et al. (2006). The angular size of NGC 2437 is well-suited to the 1 degree 2 MegaCam science 
field, and the data discussed here are the deepest wide-field images of this cluster that have 
been recorded to date. 

The (g', u' - g') and (r', g' - r') CMDs of NGC 2437 have been compared with solar 
metallicity isochrones, and a distance modulus /i = 11.05 ± 0.05 and a color excess E(B — 
V) = 0.115 ± 0.035 were obtained. These comparisons also suggest an age 220^gQ years. 
Star counts near the bright end of the r' LF favor values near the upper limit of this age 
range. The distance modulus, reddening, and age deduced for NGC 2437 from the MegaCam 
data are not significantly different from those found by Sharma et al. (2006). However, the 
distance modulus found here is larger than that computed by Stetson (1981) from multicolor 
photoelectric observations of A stars. Given its age, NGC 2437 likely formed as part of the 
recent 'local starburst' (e.g. Bonatto & Bica 2011). 

The main sequence of NGC 2437 can be traced to r' > 20 in these data, allowing the 
LF of stars with masses > 0.5 M© to be investigated. After correcting for contamination 
by non-cluster objects, the shape of the LF for stars with r' < 17 (> 0.7 M ) within 19.3 
arcmin of the cluster center is consistent with solar neighborhood star counts. However, at 
fainter magnitudes the central regions of NGC 2437 appear to be deficient in stars. 

At first glance, such a deficiency among faint stars seems to contradict the results from 
Sharma et al. (2008), who find that the mass function of NGC 2437 has an exponent that (1) 
is comparable to the Salpeter (1955) value, and (2) is consistent with that of other clusters 
in their sample. However, the faint limit of the Sharma et al. (2008) data is V ~ 18, which 
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corresponds to r' ~ 17. Their data is thus limited to the magnitude range that the MegaCam 
data finds is in agreement with the models. 

Velocity dispersion measurements suggest that there is a high binary fraction in NGC 
2437 (Kiss et al. 2008), while modelling of the cluster CMD suggests that the binary fraction 
may be 30 - 40% (Sharma et al. 2008). If stars in binary systems form in numbers that 
follow the IMF, then binarity will have the greatest impact near the faint end of the LF, 
where the incidence of binaries with mass ratios approaching unity will be highest. While the 
model LFs considered here do not include binaries, a deficiency of stars with sub-solar masses 
would still remain in NGC 2437 if models that include binaries had been used. Indeed, when 
compared with models that do not include binaries the addition of binaries typically changes 
the mass function exponent inferred from LFs by only a few tenths (e.g. Sagar & Richtler 
1991; Sharma et al. 2008). This is roughly an order of magnitude smaller than what is 
needed to explain the differences between the models and observations at faint magnitudes 
in Figure 8. 

A deficiency in stars at the faint end of the NGC 2437 LF could occur if stars with 
masses < 0.7 M Q have been dispersed from the central regions of the cluster. This could 
result from mass segregation, the structural signatures of which may be detectable after only 
a few cluster crossing times. Prisinzano et al. (2003) find evidence for mass segregation in the 
~ 100 Myr cluster NGC 2422, while Delgado et al. (2011) find evidence of mass segregation 
in even younger clusters. Sharma et al (2006) find differences in open cluster radii measured 
from optical and 2MASS data, in the sense that some clusters - including NGC 2437 - are 
more extended spatially in the near-infrared than at shorter wavelengths. Given that the 
optical and near-infrared wavelength regimes are sensitive to different stellar mass ranges, 
then this suggests that the structural properties of some clusters in their sample, including 
NGC 2437, depend on stellar mass in a way that is consistent with mass segregation. The 
investigation of the TPCF in §6 provides further evidence to support the presence of mass 
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segregation, as it indicates that the clustering properties of stars in NGC 2437 depend 
on magnitude, in the sense that bright main sequence stars have a more compact spatial 
distribution than faint main sequence stars. 

A picture is thus emerging in which NGC 2437 has a diffuse corona that contains a 
larger fraction of stars with sub-solar masses than near the cluster center. Such a diffuse 
circumcluster structure will be more susceptible to tidal stripping than the dense inner cluster 
environment, and so will likely be the source of stars that are shed from the cluster. Given 
that stars are fed into the peripheral regions of the cluster as mass segregation proceeds, 
then the corona will be replenished as stars are stripped away. In fact, the prognosis for the 
long-term survival of NGC 2437 is not good, and NGC 2437 may actively be shedding stars. 
Lamers et al. (2005) find that a cluster in the Galactic disk with a mass 10 4 M - which is 
substantially larger than NGC 2437 - may survive for only 1.3 Gyr. With a mass of ~ 10 3 
M Q (Sharma et al. 2008), then NGC 2437 will probably evaporate on the timescale of only 
a few disk crossing times. 

Is there evidence of mass loss from NGC 2437? The spatial distribution of main 
sequence stars in Figure 10 indicates that low mass stars at large radii in NGC 2437 may 
extend to the south and west of the cluster, and we suggest that this may be the roots of a 
tidal debris trail. Such trails may not be rare in the Galactic disk, and simulations suggest 
that clusters that are a few 100 Myr older than NGC 2437 may leave debris trails that extent 
over many hundreds of parsecs (Chumak & Rastorguev 2006). If the stars to the south and 
west of NGC 2437 are part of such a debris trail then they may belong to a much larger 
coherent structure, only a modest fraction of which can be captured with a single MegaCam 
pointing. 
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A. AN INVESTIGATION OF DIFFERENT CULLING CRITERIA ON THE 

CLUSTER LF 

A procedure for identifying and deleting objects with suspect photometry using e is 
described in §3, and the dataset that was culled in this manner has been adopted as the 
baseline for this study, e has the merit of being easily related to stellar magnitudes, and 
so provides - at least in principle - an intuitive basis for understanding biases that may be 
introduced by source culling. Still, e is not the only parameter that can be used to identify 
objects with suspect photometric measurements. Here, the use of two other parameters that 
are computed by ALLSTAR - sharp and x ~ to identify objects with suspect photometry 
are also explored. These parameters are described in detail by Stetson et al. (2003). They 
measure source characteristics that differ from those used to compute e, and so provide a 
means of identifying biases that may be introduced in the culling process. 

As is the case with e, sources with sharp and x that depart significantly from the mean 
value at a given magnitude may be galaxies, stellar blends, or cosmetic artifacts. Following 
the procedure employed to filter the photometric catalogue using e, a function was fit to 
the envelope of the main body of objects on the sharp vs. magnitude and x vs - magnitude 
relations. Objects that fall outside of the fitted relations were then deleted. Near the faint 
limit of the data sharp and x also show markedly greater scatter than at brighter magnitudes, 
and it was also decided to delete all objects with \sharp\ > 0.8 and x > 1-5. 

The background-subtracted r' LFs that result from cleaning the dataset using each 
parameter independently are compared in Figure Al. The faint end of the LF is shown 
in this figure, as this is where the majority of points are removed. Given that e, sharp, 
and x characterize different source properties, it is not surprising that there are LF-to- 
LF differences. However, these differences tend to be modest in size. Indeed, the LF-to- 
LF dispersion falls within the uncertainties predicted by the number counts for the LF 
constructed using e as the culling parameter. The flat nature of the r' LF at the faint end 
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is present in all three LFs, indicating that the difference between the model and the NGC 

2437 LFs at the faint end that are discussed in §5 is not an artifact of the technique used to 
excise sources with suspect photometry. 

B. THE PLANETARY NEBULA NGC 2438 

The PN NGC 2438 is located near the central regions of NGC 2437 on the sky. If NGC 

2438 is a member of NGC 2437 then this makes that cluster of considerable importance for 
studies of the evolution of intermediate mass stars and PN. However, the radial velocity 
of NGC 2438 differs significantly from that of NGC 2437, suggesting either no physical 
connection or that the PN progenitor was ejected from the cluster (Kiss et al. 2008). 

If NGC 2438 is associated with NGC 2437 then the progenitor likely had a mass that 
is comparable to those of stars in the BMS sample. When compared with lower mass stars, 
the BMS stars are tightly clustered on the sky (§6), with a typical separation of ~ 20 arcmin 
(Figure 9), which corresponds to ~ 10 parsecs at the distance of NGC 2437. Assuming virial 
eqiuilibrium, then with a total cluster mass of 10 3 M the expected velocity dispersion is 
~ 1 — 2 km sec -1 . For comparison, Kiss et al. (2008) find that the difference in velocity 
between NGC 2438 and the cluster is 30 km sec -1 . If the progenitor of NGC 2438 formed 
as a member of NGC 2437 then its high velocity (and expected large physical displacement 
from the cluster center) makes it an extremely peculiar object when compared with other 
cluster members of similar mass. 

If NGC 2438 is not associated with NGC 2437 but is part of another - not yet dissolved 
- cluster that is viewed against /through NGC 2437 then a localized concentration of stars 
may be present near NGC 2438. There are two peaks in the distribution of FMS stars in 
the lower panel of Figure 10, and the peak immediately to the north of the MegaCam field 
center coincides roughly with the location of NGC 2438. In contrast, the BMS and IMS 
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samples are more uniformly distributed near NGC 2438. Radial velocity and proper motion 
studies of objects within 5-10 arcmin of NGC 2438 (the approximate extent of the peak 
in the FMS sample) should reveal if there is an ensemble of stars that is viewed fortuitously 
through/against NGC 2437. Of course, if the progenitor of NGC 2438 is older than a few 
hundred Myr then its host cluster may have long since dissolved, although in that case the 
PN might be part of an extended debris trail. 

It is a pleasure to thank the anonymous referee for providing a timely and thoughful re- 
port that improved the manuscript. Sincere thanks are also extended to Dr. Christian Veillet 
who - as the director of the CFHT - granted the discretionary time for these observations. 
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Name 


E{B - V) 


^ log(V) 


[M/H] 


NGC 2437 


0.154 


11.17 8.4 


+0.05 



Table 1: NGC 2437 Parameters from the WEDBA Database 



a Distance modulus 
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Cluster 


Date Observed 


Exposure Times 


FWHM 


NGC 2437 


October 31, 2010 


2 sec + 95 sec (it*) 
1 sec + 50 sec (g') 
1 sec + 60 sec (r') 


1.0 arcsec 
0.9 arcsec 
0.9 arcsec 



Table 2: Log of Observations 
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Fig. 1. — The (g',u'^ — g') CMD of objects in the MegaCam science field. The envelope of 
±2(7 random errors calculated from the artificial star experiments is indicated. The open 
squares are based on the photoelectric UBV measurements of stars with V < 10.5 from 
Lynga (1960), which were transformed into the SDSS system using relations from Smith et 
al. (2002) and Jordi et al. (2006). While not the most prominent feature in the CMD, the 
main sequence of NGC 2437 runs diagonally across the CMD, and can be traced to g' ~ 20. 
The pronounced structure with u'^ — g' ~ 1.2 — 1.6 and g' > 16 is populated by main sequence 
field stars with blue colors. 



-28 - 




Fig. 2. — The same as Figure 1, but showing the (r', g' — r') CMD. The plume of sources in 
the lower right hand corner of the CMD with g' — r' — 1.5 is populated by low mass main 
sequence stars in the field. 
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12 3 

u'.-g' 

Fig. 3. — The (u^ — g',g' — r') two color diagram (TCD). UBV measurements from Lynga 
(1960), transformed into the SDSS system using relations from Smith et al. (2002) and 
Jordi et al. (2006), are plotted as open squares. The cyan line is the NGC 2437 main 
sequence traced from the CMDs, while the red line is the locus of solar neighborhood stars 
from Johnson (1966), reddened to match the mean E(B — V) of NGC 2437 and transformed 
into the SDSS system using relations from Smith et al. (2002) and Jordi et al. (2006). The 
agreement between these two sequences suggests that the u'* magnitudes track vl magnitudes 
over most of the color range. The yellow line is the ridgeline of stars in the field star plume. 
Objects with the reddest g' — r' colors at a given u'^ — g' tend to be non-cluster objects. 
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Fig. 4. — (Top Row) The (</, u'^—g') CMDs of NGC 2437 cluster members, selected according 
to location on the (u' t — g', g' — r') TCD. The three radial intervals sample equal areas on 
the sky. The main sequence of NGC 2437 is clearly evident in the inner and middle annulus 
data, but is poorly defined in the outer annulus CMD. The bluest sources in the field star 
plume have been identified erroneously as 'cluster' objects as their locations on the TCD 
fall within the extraction envelope for cluster stars. (Lower row) The CMDs of objects with 
SEDs at visible wavelengths that depart from those of cluster main sequence stars. The 
majority of non-cluster objects are faint and red. There is a more-or-less constant number 
of non-cluster objects in each radial interval, indicating that these objects are distributed 
uniformly across the field. 
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Fig. 5. — The same as Figure 4, but showing the (r',g' — r') CMDs. 
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fj. = 10.66 
E(B-V) = 0.154 
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Fig. 6. — The CMDs of cluster members in the inner annulus of NGC 2437 are compared 
with Z = 0.019 isochrones from Marigo et al. (2008). Sequences with ages 80, 160, and 320 
Myr are shown. Isochrone u' magnitudes have been transformed into u'^ magnitudes using 
the procedure described in the text. The distance modulus and reddening adopted for this 
comparison - which are those in the WEDBA database as of mid-2012 - are shown at the 
top of the figure. The models pass well above the main sequence at intermediate magnitudes, 
indicating that the apparent distance modulus is too small. 
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fj, = 11.05 
E(B-V) = 0.115 
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Fig. 7. — The same as Figure 6, but assuming an absolute distance modulus fi = 11.05 (i.e. 
r = 1.62 kpc) and E(B — V) = 0.115. Note the improved agreement between the models 
and observations when compared with Figure 6. The models are consistent with an age 
log(V) ~ 8.35 ± 0.10. 
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fj, = 11.05 
E(B-V) = 0.115 




Fig. 8. — The r' luminosity function (LF) of main sequence stars within 19.3 arcmin of the 
cluster center. N 02 is the number of main sequence stars per 0.2 magnitude interval. The 
LF has been corrected for non-cluster objects by subtracting number counts measured at 
distances in excess of 27 arcmin from the cluster center, where the density of cluster stars is 
modest (Figures 4 and 5). The error bars show uncertainties based on counting statistics. 
The green lines are solar metallicity models with a Kroupa (1998) IMF and ages 200 Myr 
(dashed line) and 300 Myr (dotted line). The models have been shifted vertically to match 
the cluster LF between r' = 13.5 and 15.5. The models (1) match the change in the slope 
of the LF near r' ~ 13, which is due to an inflexion in the mass-magnitude relation and 
provides a loose check on the adopted distance modulus, and (2) clearly exceed the observed 
star counts when r' > 17. 
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Fig. 9. — (Top panel) The TPCFs of main sequence stars in NGC 2437, selected according 
to location on the (r',g' — r') CMD. The TPCFs were constructed by dividing the separation 
functions (SFs) of objects in the BMS, IMS, and FMS samples by the SF of the FS sample, 
and then scaling the results according to the number of pairings. The TPCF of the BMS 
sample has by far the largest amplitude, and there is an inflexion point in the TPCF of the 
BMS sample at 20 - 25 arcmin separations, suggesting a change in the structure of NGC 2437 
at this point. (Lower panel) The ratio of the IMS and FMS SFs. Ratios are shown for ±0.15 
mag (solid line) and ±0.05 mag (dashed line) extraction windows along the (g' — r') axis, 
so that the effect of contamination from non-cluster sources on the TPCFs of these samples 
can be assessed. That the two ratios do not agree suggests that the relative distributions 
of the IMS and FMS samples are affected by contamination from non-cluster sources. The 
objects in the IMS and FMS samples obtained with the narrower CMD extraction window 
- which should have a lower level of non-cluster contamination - have different clustering 
properties, in a manner that is consistent with objects in the IMS sample being more centrally 
concentrated than those in the FMS sample. 
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Fig. 10. — The on-sky distribution of objects in the BMS, IMS, and FMS samples. The 
distributions of sources that fall within a ±0.15 mag extraction window about the main 
sequence locus on the (r',g' — r') CMD are shown in the top row, while the distributions 
of objects within a ±0.05 mag extraction window are shown in the bottom row. The pixel 
intensities reflect the numbers of sources in 1.5 x 1.5 arcmin 2 gathers. North is at the top, 
and East is to the left. The red circles in the middle panel of the upper row show the 
radial boundaries used in the analysis of the CMDs and the LFs. The spatial distributions 
of sources obtained with the two extraction windows are in broad agreement, suggesting 
that structural features are due to cluster stars, rather than field stars. The distribution of 
sources in the BMS sample outside of the main body of the cluster is dominated by single 
objects. The FMS sample is more diffusely distributed than the IMS sample, and has a less 
regular morphology. FMS sources extend to the south and west of the cluster, and in §7 it 
is suggested that these might be objects that have been stripped from the cluster. The PN 
NGC 2438 is located slightly above and to the left of the MegaCam field center, and there 
is a localized concentration of objects near the PN in the FMS sample that is obtained with 
the ±0.05 mag extraction window. 
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Fig. Al. — The r' LFs of main sequence stars in NGC 2437, obtained by using different 
indices computed by ALLSTAR to identify objects with suspect photometry. The error bars 
show uncertainties based on counting statistics computed in the data culled using e. The 
LFs have been corrected for non-cluster objects using the procedure described in §5. Sources 
have been excised that depart from the dominant trends with magnitude as defined by: e 
(solid line), sharp (dotted line) and x (dashed line). The LF-to-LF dispersion falls within 
the uncertainties computed for the baseline dataset discussed in §5, in which sources are 
culled using e. The flat nature of the NGC 2437 LF near the faint end is thus not due to an 
obvious bias introduced by the procedure used to remove objects with poor photometry. 



